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ABSTRACT 

Context. The high degree of deuteration observed in some prestellar cores depends on the ortho-to-para H2 ratio through the H3 
fractionation. 

Aims. We want to constrain the ortho/para H2 ratio across the LI 83 prestellar core. This is required to correctly describe the deuter- 
ation amplification phenomenon in depleted cores such as LI 83 and to relate the total (ortho+para) H^D^ abundance to the sole 
ortho-HnD* column density measurement. 

Methods. To constrain this ortho/para Ht ratio and derive its profile, we make use of the N2D+/N2H^ ratio and of the ortho- 
HiD^observations performed across the prestellar core. We use two simple chemical models limited to an almost totally depleted 
core description. New dissociative recombination and trihydrogen cation-dihydrogen reaction rates (including all isotopologues) are 
presented in this paper and included in our models. 

Results. We estimate the H2D+ ortho/para ratio in the L183 cloud, and constrain the H2 ortho/para ratio : we show that it varies across 
the prestellar core by at least an order of magnitude, being still very high (» 0.1) in most of the cloud. Our time-dependent model 
indicates that the prestellar core is presumably older than 1.5-2 x 10^ years but that it may not be much older. We also show that it 
has reached its present density only recently and that its contraction from a uniform density cloud can be constrained. 
Conclusions. A proper understanding of deuteration chemistry cannot be attained without taking into account the whole ortho/para 
family of molecular hydrogen and trihydrogen cation isotopologues as their relations are of utmost importance in the global scheme. 
Tracing the ortho/para H2 ratio should also place useful constraints on the dynamical evolution of prestellar cores. 

Key words. ISM: abundances ~ ISM: clouds - ISM: structure - Astrochemistry - Molecular processes - ISM: individual objects : 
L183 



1. Introduction 

Studies of the earliest stages of star formation have increased 
in the last ten years with the advent of new receivers acquiring 
better spatial and spectral resolution. Prestellar cores are dense 
and cold cores where gravitational collapse has not yet occured. 
In the densest regions of the core (where n^j is greater than 
a few 10"* cm"-') most heavy molecules containing carbon, 
nitrogen and oxygen seem to deplete onto the dust grains 
and only light ions remain in the gas phase. There has been 
extensive observational evidence of CO and CS depletion in the 
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center of prestellar cores (e. g . ICa selli et al 
200|, iBacmann et al] l2002l Efalla et alT 
2005) which seems to be typical of the majority of dense 
cores. Nitrogen-bearing species like CN, NH3 and NaH^ 
appear to subsist longer before freezing-out onto the dust grai ns 



?ing: 

(iTafallaet al.ll2006l iPagani et alJl2007llHilv-Blant et al.ll2008l) . 
In conditions under which heavy species are depleted, and 
(and its deuterated counterparts) are the most abundant ions 
subsisting in the gas phase. HiD^ has been wi dely detected 
and mapped in protostars and prestellar c ores (ICaselli et al.l 
I2003L IVastel etalj|2006l ICaselli et alj 1200 8') through its ortho 
fundamental line, with abundances high enough to be explained 
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by the CO depletion theory/observations. Although difficult to 
observe from Earth, the DiH^ molecule has been detected with 
its para line towa rds the 16293E pr estellar core in the L1689N 
molecular cloud jVastel et alj|2004l ) with an a bundance similar 
to the ortho-H2D^ molecule as suggested by iPhillips & Vastell 
(I2003h . 

Consequently many theoretical studies were perf ormed that in- 
cluded all the deuteral:ed form s of the ion (e.g. [Roberts et alj 
I2003L IWalmslev et all llooi) . However all nuclear spin states 
(ortho, meta, para, corresponding to the spin state of the protons 
or deuterons) must be taken into account in order to compare 
with the observational sets. Moreover the thermicity of the 
forward/backward reactions strongly depends on the symmetric 
state of the species. Though the i nfluence of the ort ho/para H2 
ratio on the chemistry of HiD^ dPagani et alJll992 j) had been 
described several years before the first detection of the ion 
jStark et al.l ll999). it is only recently that this specific sp in-state 
dependent ch emistry has been studied in detail (.Flower et alj 
I200ll2006ah . 

The motivation for our study is the many deuterated observa- 
tions performed in the L183 prestellar core (PSC) and the main 
aim of this paper is to study the ortho/para H2 ratio from some 
of these deuterated species observations. To this effect, we con- 
strain two chemical models including all the symmetric states of 
H2, D2, Hj and its deuterated counterparts with observations of 
ortho-H2D^, combined with previous N2H^ and N2D^ observa- 
tions. These models have been set up using recent dissociative 
recombination rates computed for H3 and its isotopologues as 
well as all non negligible reaction rates between H2 and H3 and 
their isotopologues (both presented for the first time in this pa- 
per). 



2. Observations 

2.1. Deuterated 

We first observed the ortho-H2D^ (and para-D2H^) with 
the Caltech Submillimeter Observatory (CSO) monopixel re- 
ceiver and subsequently took advantage of the newly in- 
stalled Heterodyne ARray Program, 16 channel 350 GHz band 
("B-band") (HARP-B) camera on the James Clerk Maxwell 
Telescope (JCMT) to fully map its emission. 



2.1.1. CSO observations 

The H2D^ and D2H^ observations were caiTied out at the CSO, 
between August 2004 and April 2005, under good weather con- 
ditions (225 GHz zenith opacity always less than 0.06). Scans 
were taken, using the chopping secondary with a throw of 
3', using the reference position: aim - 15* 54'" 08'. 50, 
<J20oo = -02° 52' 48". 

The 345 GHz (respectively 650 GHz) sidecab receiver 
with a 50 MHz acousto-optical spectrometer backend was 
used for all observations with an average velocity resolution 
of 0.04 km s"' (respectively 0.02 km s"'). At the observed 
frequencies of 372.421385(10) GHz for the H2 P+ (lio-ln) 
and 691.660483(20) for the D2H+ (Ijo-loi) lines (lAmano et al.1 
l2005i) . the CSO 10.4 meters antenna has a HPBW of about 20" 
and 11" respectively. We performed a cut in declination across 
the main dust peak in L183 between (0,-40") and (0,130") for 
H2D^ and concentrated on the reference position for D2H^. The 
system temperature was typically 1000 to 2000 K for H2D+ and 
2500 to 3500 K for D2H+. 




Fig. 1. CSO map of the H2D^ (lio-lii) line. The position is in- 
dicated in arcseconds for each spectrum and the (0,0) position 
corresponds to a2000 = 15'' 54"' 08'.50, 52000 = -02° 52' 48". 
The Y-axis represents the antenna temperature 



-2°5r00" 



-2°52'00" 



-2°53'00" 




Fig. 2. J CMT map of the H 2D+ (lio-lii) line. The dust peak 
position (iPagani et al.ll2003l) is indicated by a cross and corre- 
sponds to a20oo = 15'' 54'" 08'.50, ^2000 = - 02° 52' 48". 
Contour levels are 0. 1, 0.2 and 0.3 K km s ' 
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Table 1. Line parameters from the JCMT and CSO observations. The positions are offsets to the dust peak emission at 
aim) = 15'' 54™ 08".50, 52000 = - 02° 52' 48". For non-detected lines we give the 3cr upper limit. For JCMT and P-D2H+ at 
CSO, we give both the Monte Carlo (MC) and the LTE column density estimates. 
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" Column densities have been computed after averaging spectra at symmetrical distances from center. 

Column densities and opacities have been computed with a beam coupling correction of 70% for both JCMT and CSO. 
° Upper limit column density and opacity have been computed with a beam coupling correction of 60% at CSO. 



The beam efficiency at 372 GHz (respectively 692 GHz) was 
measured on Saturn and Jupiter and found to be ~ 60% (respec- 
tively ~ 40%) for point sources. Pointing was monitored every 
1.5 hrs and found to be better than 3". If the emission is ex- 
tended compared to the beam size of CSO, as appears to be the 
case for H2D^ then the efficiency is about 70% at 372 GHz and 
60 % at 692 GHz. The data reduction was performed using the 
CLASS program of the GILDAS software developed at IRAM 
and the Observatoire de Grenoble and the LTE data analysis with 
CASSIS developed at CESR (http://www.cassis.cesrfr). 

These CSO observations confirmed that the H2D^ emission 
was extended in this source as we could hypothesize from sev- 
eral tracers and we carried on the project at the JCMT, recently 
equipped with a new 16-pixel camera to be able to fuUy map the 
source in a reasonable amount of time. 

2.L2. JCMT observations 

The JCMT observations were obtained during semester 07A in 
service mode, using the HARP-B 16 pixel camera (one pixel, lo- 
cated in a corner, was unavailable). A third of the observations 
was obtained in jiggle-chofO mode and two thirds in position 
switch (PSw) mode. The jiggle-chop mode appeared to be no 
faster, the displacement of the telescope in PSw mode seeming 
minor compared to other overheads, and because the jiggle-chop 
mode works in the Nyquist regime, each pixel receives much less 
integration time than in PSw mode. As adjacent pixels had the 



' http://www.jach.hawaii.edu/software/jcmtot/het_obsmodes.html 



same off spectrum subtracted, the spatial average did not bring 
much improvement and we subsequently changed to position 
switch mode because deep integration on weak signal appears 
more important than Nyquist sampling for this work. In PSw 
mode, we made 2x2 pointings to fill the gaps in the camera, thus 
achieving a full beam sampling. Two such sets were performed 
to cover the main dust peak and its northern extension (Fig |2]l 
with one pixel row overlap between the two. 

Most of the observations were run in band 1 weather 
(t"225c//z < 0.05) while a few were done in band 2 weather lead- 
ing to rapid degradation of the system temperature. The source 
was observed only above 40° elevation and the band 1 weather 
system temperature was in the range 500-1000 K depending on 
the pixels and on the elevation. To observe both H2D^ and N2H^ 
(J:4-3)(at 372.672509 GHz), we tuned the receiver half way be- 
tween the two lines and used a frequency resolution of 61 kHz 
so that the backend could cover both lines at once. 

Data pre-reduction was done with the Starlink software 
(KAPPA, SMURF and STLCONVERT packages) and subse- 
quently translated into CLASS format for final reduction. 

2.2. CO depletion and dust content 

All other observations used in this paper have been obtained and 
published previously. The dust content of LI 83 both in emis- 
sion and in absorption has been described in'Pagani et alj ( I2003L 
2004). The source size is half a degree and contains a long fila- 
ment extending ~6' from north to south. Two peaks are clearly 
visible, one just south of the middle of the filament (which we 
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Table 2. Source parameters : distance from the core center, H2 
density, N2H^ abundance and N2D^/N2H^ratio. The position is 
measured aw ay from the PSC center along the R.A. axis (from 
iPagani et al.lf20(j7h . 



(60", 0) 



(45",0) 



(30", 0) 



(15",0) 



(0",0) 



Position 


H2 density 


N2H+ abundance 


N2D+/N2H+ 


(AU) (arcsec) 


'(cm-3) 









2.09 10" 


2.40E-11 


0.69 ±0.15 


1310 12 


9.23 10^ 


8.50E-11 


0.42 ± 0.05 


2620 24 


5.33 10^ 


l.lOE-10 


0.25 ± 0.02 


3930 36 


3.22 10^ 


1.53E-10 


0.16 ±0.03 


5240 48 


1.86 10= 


1.27E-10 


0.06 ± 0.02 


6550 60 


7.08 iC 


l.OOE-10 


< 0.03 



call the main peak) with an opacity of ~150 Ay and a sec- 
ond one, 3' north of the first one (the north peak) with an es- 
timated opacity of ~70 Ay. These peaks have the characteristics 
of prestellar cores. Most of the filament have an opacity above 

40 Ay. 

Two large scale C"^0 and C'^O maps obtained with the Kitt 
Peak 12-m telescope fail to trace the dense filament (Pagani et aTl 
120051) . It is now well established that this is due to depletion of 
CO onto grains. Surprisingly, the opacity at which the deple- 
tion begins (~20 Ay) is two ti mes higher than what is usually 
observed in other clouds (e.g. Alves et aLlll999l iKramer et al.l 
[1999, Bergin et al. 2002) though it still appears at a density 

1 ~3 X 10"^ fr ^*^ ^ typical depletion density threshold 

Pagani et al.ll2005l) . Possibly, the low density envelope where 
depletion has not yet occurred is very extended in this cloud 
(which is confirmed by its large influence on the C''^0 (J:2-l) 
line intensity (Pagani et al. 2 002) . The depletion f actor in vol- 
ume f or CO has been estimated to be 43 on average (iPagani et al.l 
l2005h on the line of sight of the main dust peak and is probably 
much higher in the inner part of this core where density is above 
1 X lO^* cm-\ 



2.3. N2H+ anclN2D^ 

N2 H^ and N2D^ hav e been mapped at both low (Kitt Peak 12- 
m. IPagani et al.ll2005l) and high (IRAM 30-m) resolutions. From 
the high resolution data, a strip c rossing the main d ust peak has 
been extracted and published (Paga ni et al.l 120071) . In that pa- 
per, we performed a detailed analysis of the N2H^ and N2D^ 
emission with a Monte Carlo model treating exactly the hyper- 
fine structure and line overlap of these species. We derived sev- 
eral physical properties, namely a maximum density of 2 x 10^ 
cm""*, with a radial dependence proportional to r ' up to 4000 
AU and proportional to r"^ beyond, a kinetic temperature of 7 
(+ 1) K, a slight depletion of N2H+ in the inner 3000 AU of the 
core and a deuterium fractionation which is non-measurable at 
10"* AU (< 0.03) and reaches -0.7 (± 0.15) in the center. As far 
as we know this is the highest fractionation reported yet for a 
singly deuterated species. However this may not be exceptional 
when co mpared to the detection of triply deute rated species, 
like ND. (lLiretani200l Ivan der Tak et allllool) and CD3OH 
dParise et alj|2004| ) or to the fact that most reported fractiona- 
tions are line-of-sight averages and are not derived from detailed 
profiles. 



VIsr (km/s) 



Fig. 3. ortho-H2D^ spectra across the main dust peak. East 
and west sides are folded together and fitted with a Monte 
Carlo model . Dens ity and temperature profiles are taken from 
iPagani et"an (l2007h . The spacing between spectra is 15" 



3. Analysis 

We present 3 models in this paper : a Monte Carlo model to com- 
pute H2D^ and D2H+ line intensities, a chemical steady-state 
model and a chemical time-dependent model. 

Compared to previous works, we benefit here from two new 
sets of coefficients and a large set of observations in a PSC. The 
H3 + H2 (+isotopologues) set of rate coefficients are extracted 
from the PhD work of E. Hugo in advance of publication and the 
H3 (and deuterated counterparts) recombination rates have been 
computed for this work by V. Kokoouline and C. Greene and 
are presented in Appendix IB] Rate coefficients, as computed by 
E. Hugo, describe all possible interactions between H3 and H2 
isotopologues, including reactive and non-reactive, elastic and 
inelastic collisional rates, while recombination coefficients de- 
scribe the dissociative recombination (DR) rates of trihydrogen 
cation isotopologues. Both sets take into account all the ortho, 
para and (Dp meta states. 

3. 1 . Line emission 

We have analysed the line emission of both ortho-H2D^ and 
para-D2H^ u sing a two-level Monte Carlo code (adapted from 
iBernesI 1 19791) with our new collisional coefficients. Because 
the temperature of the PSC is around 7 K for both the gas 
dPagani et al.ll2007l) and the dust dPagani et al.ll2003L |2004 . flie 
possibility of populating the next rotational level, J^y^" = 2i2 at 
1 14 K and 75 K above ground level for ortho-H2D^ and para- 
D2H^ respectively, is so low that it can be safely ignored. With 
the new coefficients, the critical densities are 1.1 x 10^ and 4.9 
X 10^ cm~^ respectively and thus the lines are close to thermal 
equilibrium in the inner core. In the case of ortho-H2D^, the 
difference between LTE (-2.0 + 0.25 x 10'^ cm"^) and Monte 
Carlo (~2.3 + 0.25 x lO'^ cm~^) column density estimates is 

typically 10-15% in the direction of the dust peak. 

The para-D2H^ line has not been detected (see lCaselli et al.l 
2008) and the 3 cr upper limit corresponds to a total column 
density of ~2.4 x lO'"' cm"^ using the Monte Carlo code. The 
LTE estimate yields -1.5 x lO'"* cm"^. 



3.2. Deuteration 

H3 ions are formed at a rate 0.96 f by cosmic ray ionization of 
H2 ( Walmslev et al. 2004), rapidly followed by reaction with an- 
other H2 to form H3 , and are destroyed in reactions with neu- 
tral species and in dissociative recombination with free elec- 
trons, negatively charged grains and possibly negatively charged 
polycyclic aromatic hydrocarbons (PAHs ). In prestellar cores, 
the primary reservoirs of hydrogen and deuteri um are H2 an d 
HD, with HD/H2 = 2(D/H)cosmic ~ 3.2 x 10'^ ( lLinskvll2007l) . 
The proton exchanging reaction of H^ with H2 partly regulates 
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HD, H„ D, H 



HCO* CO 




Fig. 4. Main reactions involved in the H3 chemical network. 
When CO and N2 are depleted, the reactions with bold arrows 
are dominant 

the H2 ortho-to-para ratio but has no effect on the abun- 
dance. Concurrently, the reaction with HD forms H2D^ and this 
primary fractionation is then followed b y the subse quent frac- 
tionations and prod uces D2H^ and tPhilhps & VasteL2003i, 
[Roberts et al]|2003l) : 



+HD 
H2D+ + HD 
D2H+ + HD 



H2D+ +H2 +232K 
D2H+ +H2 + 187 K 
D^ +H2 +234K. 



(1) 
(2) 
(3) 



The backward reactions are endothermic with an energy bar- 
rier of about 200 K (when considering only the ground level 
for each species) and were thought to be negligible at the low 
temperatures found in prestellar cores (< 20 K), in which case 



the abundance ratios 



— would be greatly enhanced. 

[''"+i^2-J„=0,I,2 

However, such enhancement can be hmited by various processes 
(see Fig.lHi : 

- dissociative recombination of H3 (and its deuterated counter- 
parts) with free electrons or negatively charged grains (and 
PAHs" ?). 

- reactions with "proton-friendly" molecules such as CO and 
N2 which destroy the trihydrogen cations to produce HCO^ 
and N2H+. 

- ortho-H2 which can react with the deuterated trihydrogen 
cation and remove the deuterium (see below). 

In this modeling, we introduce the backward reactions to 
equations (1), (2), and (3) as we distinguish between ortho, meta 
and para states of the different species. When these reactions are 
completely neglected, the deuteration fractionation is consider- 
ably enhanced and observations tow ards pre-stellar cores cannot 
be reproduced dRoberts et al.ll2003l) . Indeed, the importance of 
considering ortho and para states of various H/D carriers in the 
chemistry of the trihydrogen cat ion and isotopologues was first 
discussed by IPagani et al.l (|1992|) and subseque ntly expanded in 
a series of papers by Flower a nd coworkers (IWalmslev et al.l 
12004 iFlower et al.l2004r2006allR hereafter collectively referred 
to as FPdFW). Not only is this important in comparing the chem- 
ical model predictions on the abundance of H2D^ to the obser- 
vations of the ortho-H2D^ species alone but also because some 



important reactions are much faster with ortho-H2 than with 
para-H2 hence no longer negligible. Indeed, the 170K internal 
energy of the lowest ortho-H2 level (J=l) is large with respect 
to the temperatures of interest and can significantly enhance the 
Boltzmann factor of endothermic reactions. In some cases, re- 
actions which are endothermic with para-H2 can turn out to be 
exothermic with ortho-H2 i.e. fast and temperature independent. 
In fact, the ortho-to-para ratio of H2 is found to be a crucial 
parameter for the entire chemistry of deuterium. 

The key reactions involving ortho-H2 are essentially with 
meta-D^, para-D2H^ and ortho-H2D^ as well as para- 
H2D^(see rates in Appendix A) because the internal energy of 
the ortho-H2 alone is not enough to overcome the endother- 
micity of reactions 1 to 3. Thus only those species which also 
have an internal energy high enough (so that the sum of the 
two internal energies is higher than the endothermicity of re- 
actions 1 to 3) can react with ortho-H2 at the Langevin rate 
in cold gas. Thus, the three reactions involving ortho-H2 with 
meta-Dj, para-D2H^ and ortho-H2D^ present exothermic or 
thermoneutral channels to rehydrogenate the cations forming 
ortho-D2H^, para-H2D^ and ortho/para-H^ while the reaction 
between ortho-H2 and para-H2D^ can efficiently convert the lat- 
ter to ortho-H2D^. The ortho-H2 molecule thus opens a path 
to climb the 4 step energy ladder back from para-D2H^ to 
via para-H2D^ and ortho-H2D^ which can be very efficient in 
the presence of large ortho-H2 fractions. However, this efficient 
ladder scheme does not include D^ because the conversion of 
ortho-D2H^ to para-D2H^ is strictly forbidden in collisions with 
H2 and very inefficient in collisons with HD. Nevertheless, these 
reactions can be a strong limit to the isotopic fractionation of H3 
hence of other species. Any chemical model which includes deu- 
terium chemistry must distinguish between ortho and para states 
of dihydrogen and trihydrogen cation isotopologues and include 
reactions between the different spin states following Pagani et"aLl 
(■199Z) . , Flower et al., (,2006a.b,) and the present work. 



3.3. CO and N2 chemistry 

The CO and N2 abundances are critical parameters in the deuter- 
ation of the H3 ion. CO is expected to freeze-out onto the 
grain mantles at high de nsities (a few lO'^ c m~^) and low tem- 
peratures (< 20 K) (e . g.lCaselli e t al. 1999, Bergin et al. 200ll 
iBacmann et all l2002l iTafalla et al., .200^ .Pagani et aL .200511 
With an N2 binding energy similar to Ih e CO binding energy 
jOberg et all I2005L 'Bisscho p et all |2006|) . these two molecules 
are expected to behave similarly. However, observations towards 
prestellar cores prove that N2H^ (produced from N2) remains in 
the gas phase at higher densities than CO. This can be explained 
by the fact that N2H+ is mainly destro yed by CO (CaselH 2003, 
Pag ani et al]|2005L lAikawa et al.ll2005 ). so that the CO freeze- 
out implies a drop in the destruction rate of N2H^. This would 
partially balance the lower formation rate due to the N2 freeze- 
out. Consequently, N2H^ is observed to survive in the gas phase 



at higher densities (~10* cm ^). In the case of L183, we have 
shown that N2H^ partially survives but suffers from depletion 



at densities starting at ~5 x 10^ cm ^ to reach a factor 6^3"* at 
the core centre (~2 x 10^ cm""*). Because of growing deuterium 
fractionation, N2D+ abundance still increases towards the PSC 
center until the N2 depletion becomes predominant over the deu- 
terium enhancement, and in turn, the N2D^ abundance slightl y 
decreases in the inner most part of the core (iPagani et al.ll2007l) . 
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routes: 



The N2D^ and N2H+ ions can be produced via the following where r is the reduced temperature (t = OgrkT/e^, e being the 

electron charge, k the Boltzman constant). Therefore the recom- 
bination rate of the ion on a negatively charged grain can be 
expressed as: 



H3- 


+ N2 


^ N2H+ + H2 




(4) 


H2D+ 


+ N2 


^ N2D+ + H2 


(for 1/3) 


(5) 






^ N2H+ + HD 


(for 2/3) 


(6) 


D2H+ 


+ N2 


^ N2D+ + HD 


(for 2/3) 


(7) 






^ N2H+ + D2 


(for 1/3) 


(8) 




+ N2 


N2D+ +D2. 




(9) 



We assumed that all the H^ isotopologues react at the Langevin 
rate k^^ with N2 (which is inversely proportional to the square 
root of the reduced mass of the colliding system, hence to the 
mass of the H3 isotopologue) and that deuterium and hydrogen 
nuclei are equiprobably transferred. Consequently, H^, H2D^, 
D2H^ and D^ respectively produce an N2D^:N2H^ ratio of 0:3, 
1:2, 2:1 and 3:0. The measured ratio of 0.7 + 0.15 in the center 
of LI 83 thus imp lies significant fracti ons of D2H^ and D^ . It 
has been shown ( Walmslev et alj|2004l) that in the case of com- 
plete depletion of heavy species (C, N, O...), D3 would be the 
dominant trihydrogen cation isotopologue which would imply 
that N2D^ is more abundant than N2H^. This is not the case 
here ; nevertheless the N2H^ deuterium fractionation is a good 
constraint to the abundance of the four trihydrogen cation iso- 
topologues in our chemical model. 

At steady state (d[N2H+]/dt=0 and d[N2D+]/dt=0), reaction 
IHand its isotopic variants (|5]to|9l) being the main path to produce 
N2D+ and N2H+, we obtain : 



[N2D+] ^ [H2D+]+2[D2H+]+3[D3+] 
[N2H+] " 3[H+] +2[H2D+] + [D2H+]. 



(10) 



This ratio has been measured in the cut through the LI 83 main 
PSC. We describe in the following how our method can provide 
an estimate of the ortho/para H2D^ratio using this variable and 
subsequently, of the ortho/para H2 ratio itself as well as some 
indication of the cosmic ionization rate and mean grain size. 

3.4. Grain distribution 

Recombination of ions with electrons on negatively charged 
grain surfaces is an important process since it can be much 
faster than in the ga s phase, especially in the case of H^ 
(iDraine & SutinI [T987h . The negatively charged grain surface 
area is therefore a crucial parameter (we can safely ignore pos- 
itively and multiply negatively charged grains, consi dered to be 
very rare in cold environments. ' Draine & Sutinll987l) . The grain 
size distribution in prestellar cores is unknown since it mostly 
depe nds on grain cond ensation and also on ice condensation 
(e.g. IVastel et al.l l2006i) . We thus treat the grain radius a^v/- as 
a parameter of the model, assuming all the grains to have the 
same size and the dust to gas mass ratio to be 0.01. Different 
values could be advocated (for example, in places where deple- 
tion is importarit, the ices increase the dust mass by up to 25%, 
IWahnslev et al.l2()04h but the net result is to change only slightly 
the average grain size which is not well constrained in the PSCs 
in any case. 

The focusing effect of the Coulomb attraction between 
charged particle s and oppositely char ged grains has been in- 
cluded using the lDraine & Sutiril(ll987h formalism: 



7(Z= -1) = (1 -H - 



1 + 



2 + T 



(11) 



kgi — 



na^g,(S X 7(Z = -1)) 



(12) 



where agi is the grain radius, niu is the hydrogen mass, T the ki- 
netic temperature and S is the sticking coefficient (S < 1). The 
latter represents the probability that a colliding spec i es wil l stick 
onto the grain surface. For ions, iDraine & SutinI (|1987|) con- 
cluded that the sticking coefficient should be unity. The same 
computation can be made to estimate the recombination rate 
of other ions, H3, H2D^, D2H^, D3,- ■ ■ by a simple correction 

of the atomic mass of the ions (respectively kgj V3, koj V4, 

kg,./V5, kgr/V6, ■■•)■ In the case of collisions between charged 
particles and neutral grains, the attraction due to the polarization 
of the grain by the charged particle can be expressed through: 



/(Z = 0) = 1 -H 



(13) 



Therefore the sticking rate of electrons on neutral grains can be 
expressed as: 



k, = 



nalXS X J(Z = 0)) 



(14) 



where m^. is the electro n mass and S is t he sticking coefficient. S 
is about unity (Umebava shi & Nakanol l980) for a planar surface 
but a curvature of the grain surface will tend to reduce this pa- 
rameter However in the following we will use a factor of about 
unity as this parameter did not seem to have a large influence on 
the results in our runs. 

The grain abundance [gr] can be expressed using: 



[gr] 



(15) 



where 6 is the mean grain density (assumed to be 3 g cm"^, /d/g 
is the dust-to-gas mass ratio, and niu^ is the mass of molecu- 
lar hydrogen. Another important parameter in our model is the 
abundance ofthe negatively charged grains ([gr] = [gr°] + [gr"]). 
At steady-state, assuming partial depletion of CO and N2 and to- 
tal depletion of all the other heavy species : 

d[gr-] 



dt 



[gr"][e-]^e - [gnimksr - [gr-][Hj]^g,3 

- [gr-][H2D+]^g,4 - [gr-][D2H+]^g,5 - [gr-][D3+]^g,6 

- [gr-][N2H+]^N.m - [gr-][N2D+]A:N,D^ 

- [gr-][HCO+]^HCO- - [gr-][DCO+]^DCO- = 

(here we have neglected HD^, Dj, He^,...) 



3.5. Steady-state cliemical model 

The code we describe in the following is used to calculate the 
steady-state abundances of the chemical species found in the dif- 
ferent layers of the LI 83 prestellar core as listed in Table |2] 

In the steady-state approximation the abundance species are 
interlinked via their production rates and their destruction rates 
(production=destruction). 

Since Ht is produced at a rate 0.96^, the H^ abundance can be 
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expressed as (including only the main reactions) : 

0.04^ 



[H-] 



(16) 



The main production path of Hj is via cosmic ray ionization 
of H2 and proceeds in two steps : 



0-H2 + e 



^ + 0-H2 
( + P-H2 P-H2 + e 



(17) 
(18) 



and rapidly reacts with another H2 to form H3 but the branch- 
ing ratios between different combinations of spin states are non- 
trivial dOka 2004) : 



P-H+ + P-H2 ^ P-H+ + H 

p-H^ + 0-H2 ^ P-H+ -H H 2/3 

^o-HJ-hH 1/3 

o-Hj + P-H2 ^ p-Hj +H 2/3 

^ o-Hj -H H 1/3 

0-H+ + 0-H2 ^ P-H+ -H H 1/3 

^o-H+-hH 2/3 



(19) 
(20) 

(21) 

(22) 



These are different from those advocated by FPdFW who took 
branching ratios of 1/2 for both species. The ortho-H^ forma- 
tion rate from cosmic ray ionization k^r-o is therefore the sum of 
several terms : 



0.96(l/3[p-H2][o-Hj] 



l/3[o-H2][p-Hj] + 
2/3[o-H2][o-Hj]) 



(23) 



The production rate for ortho-H^ can be expressed as (including 
only the main reactions. The rates are listed in Table lA.ll : 

kcv^^ + ^-loood[o-H2][o-H2D^]nH2 + (^Opoop + 

^0pooo)[p-H+][o-H2]«H, + ^lpdod[HD][p-H3+]nH, (24) 

which represent respectively the formation from cosmic ray 
ionization, backward destruction of ortho-H2D^with ortho-H2, 
spin conversion of para-Hj with ortho-H2 and finally, spin con- 
version of para-Hj with HD. 

The destruction rate for ortho-H^ can be expressed as (in- 
cluding only the main reactions) : 

(0-^recl[e"] + (^Ooopp + ^Ooopo) [0-H2] 
+ (fclodpd + ^lodpo + ^lodop + ^lodoo)[HD] 

+k,,[CO] + kmim + kg,i[gr-]))[o-H+]nH, 



(25) 



which respectively represents its destruction by dissociative re- 
combination with electrons, spin conversion with ortho-H2, spin 
conversion and deuteration with HD, proton transfer reactions 
with CO and N2 and dissociative recombination on grains. 

Similarly, the para-H^ formation from cosmic ray ionization 
can be expressed as : 



^„-p = 0.96([p-H2][p-H^] + 2/3[p-H2][o-H^] + 
2/3[o-H2][p-H+] + l/3[o-H2][o-H+]) 



(28) 



The production rate for para-H^ is : 

kcr-p( + (A:_loopd[o-H2][0-H2D+] + 

(A:Ooopp + ^0oopo)[o-Hj][o-H2] + A:lodpd[HD][o-H+] + 

^_lpopd[o-H2][p-H2D+])«H, (27) 
and the destruction rate is : 

(f-^recl[e"] + (^Opoop + ^0pooo)[o-H2] 
+(^lpdod + ^Ipdpo + ^Ipdop + ^lpdoo)[HD] 

-H^eo[CO] + kmim + ^gii[gi-"]))[p-Hj]«H2 

The N2 abundance has been solved numerically to obtain the 
observed N2H^ abundance. Electronic abundance is adjusted to 
reach equilibrium. 

In our steady-state model, the H2 ortho/para ratio, the aver- 
age grain radius and the cosmic ionization rate ( are the varying 
input parameters. Within each layer of the PSC model (Table|2]l, 
these parameters are adjusted to match the following : 

- the H2 density 

- the N2D+/N2H+ ratio at 7K 

- the observed ortho-H2D^ column density 

- the upper limit on the P-D2H+ column density 

Although the full range of grain sizes and ortho-to-para H2 ra- 
tios have been explored for each H2 density, we have not al- 
lowed solutions in which, for example, the grain size would os- 
cillate from one layer to the next. We have searched for solutions 
throughout the layers in two different scenarios : (i) the grain size 
and the ortho-to-para ratio were both kept constant ; or (ii) the 
grain size was allowed to increase and the ortho-to-para H2 ra- 
tio to decrease with the H2 density. In both scenarios, ^ was kept 
constant througho ut the l ayers. We neglected detailed reactions 
with D2 as Flower et al.l ( |2004|) have shown that its role is neg- 
ligible in general and we have kept the HD abundance constant 
which is generally a good approximation. 



3.6. Time-dependent chemistry 

In a second step, we have constructed a pseudo time-dependent 
model based on NAHOON, a chemical model, a version of 
which has been made publicly available by V. WakelarrQ- We 
have modified this model in two ways : 1) we have replaced 
electron (resp. ion) reactions with neutral (resp. charged) grains 
as provided in the Ohio State University (OSU) reaction file (de- 
livered with NAHOON) by the set of equations described above 
(Sect. [34l and l3.5l l, which we have directly included in the pro- 
gram, to take into account Coulomb focusing ; 2) we have in- 
cluded the formation of HD and D2 on grain surfaces and we 
have introduced the spin state of H2 and D2 with the usual as- 
sumption that they are formed with the statistical ortho/para spin 
state ratio of 3 and 2 respectively. We have used the formation 
rate provided in the OSU reaction file (5x 10 cm^s ') for 
the formation of molecular hydrogen. Because grains are cov- 
ered by ice in the environments concerned here, we consider 
that the only interaction between the atoms and the surface is 
physisorption. In this case, the formation rates of HD and D2 (in 
cm^ s ') is half for HD and 10^ t imes lower for D2 w ith respect 
to H2 formation (as dicussed in iLipshtat et al.l 12004 ). In envi- 
ronments where grains are not covered by icy mantles, on the 
other hand, one would have to consider chemisorption, which 



(26) 



- http://www.obs. u-bordeaux 1 .fr/radio/VWakelani/Valentine%20 
Wakelam/Downloads .html 
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strongly changes the efficiencies of the formation of HD and D2 
dCazaux et al.ll2008h We have reduced the set of species and re- 
actions to our needs, limiting ourselves to the most important 
reactions (see below) but differentiating all ortho and para (plus 
meta-Dp species as independent species and including all the 
detailed rates between the trihydrogen cation and dihydrogen 
isotopologues (see Sect. I3.7.T1 ) including spin state conversions. 
However we have included more reactions than in the steady- 
state model, taking into account reactions with D2, , He^, etc. 
and allowing the ortho/para H2 ratio and the HD abundance to 
vary. 

The main path to convert ortho-H2 into para-H2 is via the 
reaction 



0-H2 -H H+ ^ P-H2 + H"^ 



(29) 



which proceeds seven orders of magnitude faster at 7 K than the 
reverse reaction. 



3. 7. Rate coefficients 

Many groups have made available gas-phase rate coeffi- 
cients. The University of Manchester Institute of Science and 
Technology (UMIST) Database for Astrochemistry contains in- 
formation on 4500 reactions of which 35% have been mea- 
sured experirnentally , some at temperatures down to 20 K 
(IWoodall et al.ll2007h . The OSU group provides approximately 
the same database but focuses more on low temperature chem- 
istry. We accordingly use in our modeling some of the reactions 
in the latter (with the most recent version OSU2007), consider- 
ing the low temperatures found in L183. Apart from the rates 
presented below, all reaction rates involving deuterium have 
been taken from FPdFW (except recombination on grains for 
which we have used different sticking probabilities). We disre- 
gard the odd branching rat i o of N ^H^ dissociative recombination 
reported by Geppert et al.l (|2004|) to consider a single p ossibility, 
namely the liberation of dinitrogen dMolek et al.ll2007h . 



3.7.1. H3 + H2 isotopologue reaction rates 

Phase space theory (PST) was used to derive thermal state-to- 
state rate coefficients for the whole H^ + H2 — > H^ -H H2 sys- 
tem and isotopic variants in the temperature range 5-50K. This 
statistical method accounts for such quantities as mass, energy, 
rotational angular momentum, nuclear spin symmetry and their 
respective conservation laws. The ergodic hypothesis which is a 
requisite for PST as well as the full-scrambling hypoth esis are 
assumed according to the topology of the PES ( Yama guchi et al.l 
ll987llXie et al .l2005h . Reactant (product) trajectories are treated 
with the classical Langevin model. The resulting set of state-to- 
state rate coefficients deviates from the detailed balance princi- 
ple by a few percent at worst and is consistent with thermody- 
namical equilibrium constants. Details will be given in a forth- 
coming publication (Hugo et al., in prep.) 

In the present astrochemical model, nuclear spin states of the 
different molecules are treated as distinct species but their rota- 
tional states are not considered individually. We thus made the 
assumption that only the rotational ground states of each nuclear 
spin species were populated and used the ground state-to-species 
thermal rate coefficients obtained by summing ground state-to- 
state thermal rate coefficients over the product channels. 



3.7.2. H3 isotopologue dissociative recombination rates 

The dissociative recombination (DR) rate of H3 is diffi- 
cult to meas ure exper imentally. Since the early measurements 
e.g. by A dams et al.l {1984), numerous att empts have been 
made and are summarized in two papers jBates et al] I1993L 
iFonseca dos Santos et al.l l2007h . In parallel, theoretical work 
has also been develope d with the latest one published by 
Fonsec a dos Santos et al. I (Hoot"). Extending upon that work, 
we present here, in Tab le IB. 11 the H^ updated DR rate 
dFonseca dos Santos et al.l 120071) along with newly calculated 
H2D^, D2H^, and D^ DR rates (see appendix B for more de- 
tails). These calculations do not predict the branching ratio of 
the DR products. We have thus adopted the branching ratios 
jjublished elsewhere CDatz et al . 1995a b, Strasser et al. 200j, 
Zh aunerchvk et al. which we have applied to the calcu- 

lated rates. The resulting DR rates at 7 K are listed in Table lA.fl 
Several remarkable effects at low temperature are visible (see 

Figs 10 121 13 : 

- the strong departure of the ortho-H^ DR rate from that of 
para-Hg . The difference is a factor of 10 at 10 K. 

- D2H^ shows a large DR rate drop, by a factor of 10 at 10 K 
for both ortho and para species compared to the extrapolated 
value used by FPdFW. 

- On the contrary, a large increase of the D^ DR rate is pre- 
dicted to occur but mostly at temperatures where deuteration 
is low and therefore the consequences for D3 abundance are 
limited. 



4. Results and discussion 

4.1. H2D+ linewidtii 

In order to fit the observed ortho-H2D^ line profiles, we run the 
Monte Carlo model using the "best model" velocity, density and 
tempera ture profiles derived from the N2H^ and N2D^ data ana- 
lyzed in lPagani et al.l d2007h . However, the linewidth for the two 
stronger spectra (offsets (0,0) and (15", 0), Fig.|3]l is too wide to 
be reproduced with the same micro-turbulent width which we 
have used for N2H+(Avtuib(FWHM) ^ O.Ukms '). Indeed, the 
central spectrum linewidth measured by fitting a Gaussian yields 
0.45 (+0.03) km s"'. The thermal contribution is 

'kT 



Av,herm(FWHM) = 2.336 X J — = 0.28 km s"' (30) 

V m 

at 7 K, which implies that Avtu,b(FWHM) contribution is 0.35 
km s 2.5 times larger than for N2H^. Here k is the Boltzmann 
constant and m is the mass of H2D^. If we impose a turbulent 
velocity similar to the one modeled for N2H^, then the temper- 
ature needed to obtain such a wide line wo uld be 16K which i s 
completely ruled out by N2H^ observations (Pa gani et al.ll2007l) . 
Infall motion limited to the inner core could have explained the 
H2D^ width if no H2D^ spectrum other than the central one had 
been wide combined with a large depletion of N2H^. This is not 
possible because the H2D^ spectrum at (15", 0) has the largest 
width (0.49 +0.03 km s"') in a region where N2H^ is hardly 
depleted. This remains therefore a pending problem. 



4.2. A/2H+ deuteration 
4.2.1. Requested conditions 

We next discuss the main parameters that control the N2H^ 
deuteration using the steady-state chemical model. 
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Fig. 5. N2D^/N2H^ ratio as a function of the ortho/para H2 ratio for all possible electronic abundances and total depletion (the 
steady-state chemical model). The lower row corresponds to the densest part of the PSC (n(H2) = 2 x 10^ cm"^) and the two 
horizontal dashed lines the measured range of the N2D^/N2H^ ratio while the upper row corresponds to the external part of the PSC 
(n(H2) = 7x10"^ cm"^) with the dashed line representing the N2D^/N2H^ ratio upper limit. The three columns represent different ^ 
values as indicated above {(0 - 1 x 10"'^ s"') . The large arrow indicates the direction of increasing electronic abundance 



? = 10 X 



f = 100 X 



10" 

10^ 

q 

E 10^ 

+ 
z 

*Q 10"^ 

-z. 

10-^ 
10-^ 
10^ 



10 



10 



Q 10 



10 



10 



!!!!)!, 



[el 



n(H2) = 7.0 X 10* cm ' " 



-|— I — h- 1— I — h- 1— I — h- 1— I — h- 1— I — h- 1— I — 1-4- 



1 n ! n 



[el 



n(H2) = 2.0 X 10^ cm ^ 
. . . . I . . . . . I 



\ \ \ \ 



[el 



H h- 1— I Y 



-l— I — h- 1— I — h- 1— I — h- 1— I — h- 1- 



[el 



5 \ 



[el 



-i— I — h- 1— I — h- 1— I — h- 1— I — 1 1 1 



[el 



10 



-6 



10 



10 



1 10 



10 10 
H2 ortho/para ratio 



1 10 



10 



10 



Fig. 6. Same as Fig. |5]but with a CO/H2 abundance of 10 ^ in the outer layer (n(H2) = 7 x 10"^ cm ^) and 10 ^ in the inner layer 
(n(H2) = 2 X 10*^cm-3) 



The models have been run for a temperature of 7 K which 
prevails in all the laye rs where N2D'^ has been detected in the 
PSC cut presented in 'Pag ani et al.l (|2007|) . We have also run 
the models for the corresponding density, N2H+ abundance and 
N2D^/N2H^ ratio of each layer (the parameters are listed in 
Table©. 



As discussed above, the abundance of ortho-H2 is the main 
controlling factor of the trihydrogen cation isotopologue abun- 
dances and therefore of the N2 D'^/N2H^ ra tio (and similarly of 
the DCO+/HCO+ ratio, see e.g. lMaret &"B ergin 2007). We have 
therefore explored the range of possible solutions for the or- 
tho/para H2 ratio in the two extreme layers of our core profile 
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(n(H2) -1 y. 10"* cm ^ and 2 x lO^cm ^) for which we have a 
N2DVN2H+ ratio of <0.03 and 0.7 + 0.15 respectively (Fig.|5l). 
We have done this for three cosmic ray ionization rates (10"'^, 
10 and 10 s~^ ) covering the values generally discussed 
in the literature (e.g. lMaret & Berginll2007l) and for all possible 
electronic abundances (or average grain sizes as they are linked 
via the abundance of which is mostly controlled by the grain 
surface area). In this first run, we have simulated total depletion 
by adjusting the CO and N2 abundance|3 to 10"^. We have also 
indicated the range of N2D^/N2H^ ratio measured in both lay- 
ers. The average grain radius has been varied from 0.01 /vm to 
5 and electronic abundance from 10"" to 10"^ which cover 
the usually accepted values. We can see that N2D^/N2H^ ratios 
above 100 are possible in dense gas though they require very low 
electronic abundances and therefore very small grains which are 
probably absent fro m these dense and co ld regions due to grain 
coagulation (see e.g. lStepnik et al.ir2003h . 

In the lower density outer layer where no N2D^ has been de- 
tected, the ortho/para H2 ratio must be high enough, i.e. above 
~0.01, to prevent any deuteration from occuring whatever the 
cosmic ray ionization rate. On the contrary, the dense, strongly 
deuterated layer has solutions only below a maximum ortho/para 
H2 ratio of 0.01 (or lower for high ^ rates). Thus the ortho/para 
H2 ratio across the PSC clearly must vary from above 0.01 to 
below 0.01. In the case of low cosmic ray ionization rates (10"'^ 
s"'), though the ortho/para H2 ratio of 0.01 seems to be a com- 
mon solution for both layers, it requires a high electronic abun- 
dance (and large grains) in the outer layer and a low electronic 
abundance (and small grains) in the inner dense part. This is 
clearly improbable. The temperature being low enough in all the 
layers, warm layers (above 20 K) cannot be invoked instead of 
a high ortho/para H2 ratio to limit the deuteration in the outer 
parts of the PSC. CO total depletion is however questionable and 
in the model we also used a CO depletion factor of 10 (abun- 
dance of 10"^) in the outer layer and a CO depletion factor of 
100 (abundance of 10"^) in the inner layer (Fig. |6]l. This only 
limits the maximum N2D^/N2H^ ratio which decreases by one 
order of magnitude. Indeed, the destruction of by CO domi- 
nates over recombination with electrons when their abundance 
is very low and vice versa. However, the conditions to reach 
the observed N2D^/N2H^ ratio remain unchanged and therefore 
only a variable ortho/para H2 ratio can be invoked. Such a vari- 
able ortho/para H2 ratio cannot be investigated with a steady- 
state model because in all layers, the ortho-H2 abundance would 
eventually decrease to values about of lO^-'-lO"'* as discussed by 
FPdFW. 

4.2.2. The ortho/para H2 variation 

We discuss here the possibilities of making the ortho/para H2 
ratio vary across a single PSC. 

It is commonly accepted that H2 is formed on grain surfaces 
with an ortho/para ratio of 3 because of spin statistics and the 
exothermicity of the reaction H + H — > H2. Subsequently, the 
ortho-H2 is converted into para-H2 following equation|29]and to 
a lesser extent with reactions involving and its isotopologues. 
As already discussed by FPdFW, this conversion is slow and has 
probably not reached steady state in clouds with ages between 
10^ and 10* years. 



' Here, we look for general solutions of the N2D^/N2H^ ratio as a 
function of several parameters, we thus do not try to fit the N2 abun- 
dance to obtain the observed N2H+ abundance 



We have therefore run the modified Nahoon model to search 
for a simultaneous solution for all layers which would get the 
necessary ortho/para H2 gradient across the PSC to achieve both 
the observed N2D^/N2H^ ratio and ortho-H2D^ abundance pro- 
files. In this model, we have no direct measurement of the CO 
abundance but for the credibility of the model, we have set the 
CO abundance to 10"^ in the outer layer and increased the deple- 
tion to reach a factor of 100 (i.e. a CO abundance of 10"*) in the 
core. The results are presented in Figs.|7]&[8] We have searched 
for a solution where each modeled layer meets the two obser- 
vational constraints (N2D^/N2H^ ratio and ortho-H2D^ abun- 
dance) at the same time, but these solutions are not simultaneous 
between the different layers. We varied the average grain radius 
and the cosmic ray ionization rate, ^. We could find solutions 
for grains of average radius 0.025 to 0.3 fim . No solution has 
been found for grains above 0.3 yum. Figs. |7]&|8] show the case 
for which the grain average radius is 0.1 fim and f = 2 x 10"'^ 
s"'. In this case, the time range inside which all layers meet the 
requested conditions is 0.6 to 1.7 x 10^ years. Figure [8] shows 
how the ortho/para H2 ratio evolves for 3 selected layers. We 
have marked the appropriate time which is the solution for each 
of these layers as established from Fig. Q In that figure, we can 
see that the o/p H2 ratio is below 0.01 for the dense layers and 
still above 0.01 for the outer layer for which no H2D^ has been 
detected, as expected from the steady-state model. We can also 
see that the full o/p H2 relaxation has not yet occurred even for 
the densest part of the cloud. Smaller grains have a larger in- 
teracting surface and therefore lower the abundance of ions 
which preferentially recombine on negatively charged grains (or 
PAHs"), consequently slowing down the dominant ortho-H2 re- 
laxation reaction (eq. |29] l. Though smaller grains also imply a 
lower electronic density, therefore favouring a higher deutera- 
tion of N2H^ as shown in Fig. |5] the slower disappearance of 
ortho-H2 is the dominant process here and finally, smaller grains 
slow down the deuteration process. For grains of average radius 
0.025 /im, the range of ages matching the range of N2H^/N2D^ 
observed ratios is 2.7-3.8 x 10^ years, while for grains of aver- 
age radius 0.3 fim, the time range is only 3.8-7.2 x lO'* years. 
Figure |8] suggests that D2 should become a sizeable fraction of 
available deuterium a short while after the present state (typi- 
cally 2-3 X 10^ years) and that HD coiTespondingly should drop 
slightly. 

4.3. Age of the core and collapse 

Though it is normal that dense layers evolve faster than less 
dense ones, at least to account for a differential ortho/para H2 
ratio, the layers are evolving too fast in our model. The dens- 
est layers would have reached their present status 2 to 3 times 
faster than the outer, less dense layers. This could be possible 
only if the denser layers had reached their steady-state equilib- 
rium. This is not the case here where the densest layer would 
reach its steady-state equilibrium only after 2 x 10^ years and 
this would imply a N2D^/N2H^ ratio of 6, almost an order of 
magnitude higher than observed. The most probable reason for 
this time discrepancy is that the core has undergone a contrac- 
tion and therefore all layers were not so dense in the past. While 
the outer layers have only little evolved in density (the most ex- 
ternal one probably started at 0.5-1 x 10^ cm"-' to reach 7 x 10^ 
cm"^ today), the inner ones have undergone a much greater den- 
sity increase. As constant density through the core would give 
no chemical differentiation while a time-frozen density profile 
as measured here gives too much differentiation, the solution is 
in between the two. Starting from a uniform gas, the chemical 
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differential evolution of the core should therefore help us to con- 
strain the duration of the contraction and the type of contraction. 
Of course the model should also include the evolution of deple- 
tion which also plays a role in the acceleration of the process. 



As the core must have started to contract from a lower den- 
sity region, typically 10"* cm"^, it is clear that all layers have 
accelerated their chemical evolution while their density was in- 
creasing. Therefore, the layer with the longest time to reach the 
observed N2D^/N2H^ ratio gives a lower limit to the age of the 
cloud. Depending on the exact average size of the grains, this 
lower limit is 1.5-2 x 10^ years here. It is even larger because 
before the cloud underwent contraction, depletion had not yet 
occured and therefore species like atomic sulfur, S, must have 
been present in quantities large enough to transfer notable quan- 
tities of electric charges from to (H^ -i- S — > H -1- S^) and 
PAHs must have also be en abundant enough to help destroy 
ions jWakelam & Herbstil2008l) . All these phenomena contribute 
to the diminution of the H^abundance, th erefore slowing down 
the ortho-H2 relaxation process. Indeed. [Flower et al] ( 2006bh 
show that the relaxation process in some cases takes 3 x 10^ 
years, typically 50 times slower than in the case presented here 
(and 15 times slower for similar conditions of grain size and cos- 
mic ray ionization rate but without depletion in their case). 



4.4. Para-DzH^ 

At 7 K, the strongest possible line intensity (LTE case) for the 
ground transition of para-D2H^ is below 0.3 K because of the 
Rayleigh- Jeans correction at 691 GHz which becomes large. 
Moreover the thermalization of the line is difficult to obtain in 
this source beyond the radius of 3000 A.U. because the PSC 
density drops below the para-D2H^critical density (nc„, = 4.9 
X 10^ cm"^) and a slight drop of the excitation temperature turns 
into an exponential decrease of the brightness temperature. For 
Tex = 6 K, T bright < 0.13 K. Therefore searches for para-D2H^ 
must reach very low noise levels to have a chance of detection. 
From the models we present here, we predict an integrated line 
intensity of 1 1 mK km/s (34 mK peak) with the Monte Carlo 
model, and an upper limit of 16 mK km/s (48 mK peak) in the 
case of LTE. This is a factor of 3 to 4 below the upper limit we 
obtained from the observations. 



4.5. The chemical profile 

Finally, we obtain a detailed profile of the PSC which we present 
in Fig. |9] It represents the solution for the model we presented 
in Figs. |7] & ID taking for each layer the values at their respec- 
tive best-fit time. The large variation of the ortho-H2 species 
across the core (a factor of 15) makes change by an even 
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larger amount (2 orders of magnitude) but it does not become 
the most abundant trihydrogen cation isotopologue in the core 
center at this stage because the ortho-H2 abundance is not yet 
low enough. In the present case, for a density of 2 x 10^ cm"^, 
the inversion between H3 and occurs when the abundance of 
ortho-H2 drops below 3 x 10"^. 

The ortho/para H2 ratio in the outer layer is 0.04 (as expected 
from Fig. |6] which indicates a lower limit of 0.01). As discussed 
above, the ortho/para H2 ratio evolution speed is linked to den- 
sity and grain size, both of which are lower outside the PSC, in 
its embedding parental cloud. We can thus expect this ratio to be 
at least 0.05 and probably above 0. 1 in the envelope of the cloud. 

N2 is an input parameter in our model because we have 
not included all the nitrogen chemistry. As reactions with 

isotopologues ar e the main path to destroy this molecule 
jFlower et al." 2006b), we do not introduce a large error in deter- 
mining its abundance directly from the N2H^ abundance itself 
and obtain a N2 profile which is probably closer to reaUty than 
if we had let the whole N-chemistry freely establish its abun- 
dance because of too many unknowns. The abundance profile 
thus starts at 1.5 x 10"^ with respect to H2 in the low density 
layer to diminish to 3 x 10"^ in the densest layer. The unde- 
pleted N2 abundance after attaining the steady state is ~3 x 10"^ 
( iFlower et all 12006 b') but this is reached only after » 5 x 10'' 
years. As depletion of N2 in the outer layer is possibly still small, 
we can conclude that N2 has not yet reached its steady state 
abundance which puts an upper limit on the age of the cloud 



of abou t 1 X 10^ years following the estimate of iFlower et al.l 
(l2006bl) . However, this depends very much upon several factors, 
for example the C:0 elemental abundance ratio. Consequently 
this information is only indicative. 

Though our model does not deal with the nitrogen chemistry, 
it seems to indicate that low abundances of N2 are sufficient to 
explain the observed N2H^ abundance and therefore, the N2 de- 
pletion seems not to be a critical factor as long as CO is also de- 
pleted. The much debated contradiction between the presence of 
N2H^ in depleted cores while N2 should deplete like CO would 
thus be a false problem. Low CO and electronic abundances, 
limiting the destruction rate of N2H^, seem to be sufficient to 
compensate for the N2 depletion itself to a large extent. Finally, 
in the inner core, N2 and N2H^ (+ N2D^) follow a similar de- 
creasing trend, suggesting that N2 depletion eventually forces 
N2H^ to decrease. 

While ortho-H2D^ is 83 K above the para ground state, it 
is more abundant all over the core profile by almost an order 
of magnitude whereas the thermal equilibrium ratio would be 
ortho/para H2D^ a; 2 x 10"^ at 7 K. This demonstrates the ef- 
ficiency of the ortho-H2 in converting para-H2D^ into ortho- 
H2D^and therefore limits the total abundance of H2D^, as the 
backward channel to H3 remains open even at 7 K. This or- 
tho/para population inversion does not occur for D2H^ as the 
species needed to perform this inversion is no longer ortho-H2 
but the much rarer ortho-D2. Therefore, the para-D2H^ remains 
the least abundant of the two spin state species which, combined 
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with the fact that its ground transition is higher in frequency than 
the one of ortho-HaD^, makes its detection extremely difficult. 

5. Conclusions 

We have presented a pair of simple chemical models restricted 
to H-carriers, He plus CO and N2 to account for the observed 
HCO+, DCO+ (not discussed in this paper), N2H+, and N2D+ 
ions. We have benefited from new computed reactions rates for 
both the + H2 isotopologue combinations and for the 
isotopologue dissociative recombination rates which explicitly 
take into account the nuclear states individually. 

With the steady-state model we have shown that the or- 
tho/para ratio of H2 must vary from above 0.01 in the outer parts 
of the L183 PSC to less than 0.01 in the inner parts to explain 
the variation of deuteration across the core. Checking the reality 
of the ortho/para H2 variation with a time-dependent model , we 
have also found that if the present PSC density profile is static, 
then the inner layer would have reached its present status 2 to 3 
times faster than the outer layers. Because the present status is 
not in steady-state, the layers should evolve at a similar rate and 
therefore the density must have been lower in the past. The most 
probable explanation is that the core has probably evolved from 
a uniform density cloud to the present centrally condensed PSC. 
The time-dependent model also suggests that the ortho/para H2 
ratio changes by one order of magnitude from ~5% at a density 
of 7 X 10"* cm"^ down to a few x 10"^ in the inner dense core. 
This has two important consequences : 

- it is most probable that most of the cloud, outside the densest 
regions (i.e. the two PSC and the ridge in between) have an 
ortho/para H2 ratio also above 5%, and possibly 10%, con- 
trarily to what is usually assumed in models. 

- In principle, it should be possible to fit the PSC profile with 
this chemical model combined with a dynamical model in- 
cluding depletion, to set an age to this PSC and possibly dis- 
criminate between several types of collapse but it is beyond 
the scope of this paper. 

We already have some indication that the age of the PSC is 
somewhat above 1.5-2 x 10^ years though the N2 abundance 
suggests a relatively short time (10^ years, except if depletion is 
compensating for its formation) when we compare our adjusted 
abundance to the formation rate of N2 given by Flower et alj 
(l200 6b). The low abundance of N2 needed to explain the ob- 
served N2H^ abundance indicates that its depletion is not a real 
problem, though, obviously, N2H^ would be much more abun- 
dant if N2 was not being depleted but CO still was. 

Finally, we stress the importance of considering or- 
tho/meta/para chemistry when dealing with the deuteration of 
the interstellar medium. The importance of the ortho-H2 on 
the amount of deuteration and the observations limited to only 
ortho-H2D^ make this inclusion compulsory. Moreover, a com- 
plete state-to-state chemical model should be developed to take 
into account rotational pumping, leading to a higher destruction 
rate of the deuterated trihydrogen cation and possibly explaining 
the observed linewidth of ortho-H2D^. 

Detecting para-D2H^ would be highly desirable to help con- 
strain the models, but the high frequency and limited trans- 
parency of the atmosphere make it a difficult tool to use. Though 
observable from the ground, because of its weakness in cold dark 
clouds, which are the only places where it should be found, di- 
rect para-D2H^ observations should be made on a large number 
of Galactic lines of sight using the HIFI receiver on board the 
Heschel Space Observatory. 
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Appendix A: The reaction rate table used in the Nahoon modified chemical model 
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Reaction 



Rate 



Reference 



Table A.l. Reactions used in the model. The rate coefficients are given for 7K. 



Reference (1) corresponds to lGerlictil il990h . references (2), (3), and (4) corre- 
spond to E. Hugo, OSU 07, and this paper respectively. For OSU 07, branching 
ratios involving spin states have been infered from quantum mechanical rules. 
For reactions involving grains, a gr ain radius of 0.1 ii m a nd a sticking coeffi - 
cient of 1 have been considered. (5) 'Datz et al 



(7) IZhaunerchvk et al.l j2008l) . (8)|^sson et al 



1995h^ (6 }lDatzetalJ ( ll995j) . 
19971). (9) Riolek et alj]200l . 
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Appendix B: On the rate coefficients for 
dissociation recombination of IH^ and its 
isotopologues 

The dissociative recombination (DR) r ate coefficients for ortho- 
and p ara-H^ have been published by iFonseca dos Santos et alj 
( 12007 ). Here, we present the results obtained for all four iso- 
topologues. The DR rate coefficients for different species of the 
nuclear spin are calculated using the approach described in a se- 
ries of papers devoted to DR theory for triatomic molecular ions. 
See iKokoouUne & Greend (2003a b), Fonseca dos Santos et al. 
(2007) for H3 and D^ calculations and Kok oouline & Greene! 
(2004, 2005) for H2D+ and D2H+. The scope of this paper does 
not allow us to review ffie theoretical approach in detail. We only 
list its main ingredients. 

The theoretical approach is fully quantum mechanical and 
incorporates no adjustable parameters. It relies on ab initio cal- 
culations of potential surfaces for the ground electronic state of 
the H3 ion an d several excited stat es of the neutral molecule H3, 
performed bv lMistrik et al.l (1200 ll) . 

The total wave function of the system is constructed by 
an appropriate symmetrization of products of vibrational, rota- 
tional, electronic, and nuclear spin factors. Therefore, rovibronic 
and nuclear spin degrees of freedom are explicitly taken into ac- 
count. 

The electronic Born-Oppenheimer potentials for ffie four H3 
(and H3) isotopologues have the C3,. symmetry group. The C3,i 
symmetry group has a two-dimensional irreducible representa- 
tion E. The ion has a closed electronic shell. The lowest elec- 
tronic state of the outer electron in H3 has p-wave character. The 
;:i-wave state of the electron also bel ongs to ffie E representa - 
tion. Due to the Jahn-Teller theorem dLandau & Lifshit3l2003l) . 
this leads to a strong non-adiabatic coupling between the E- 
degenerate vibrational modes of the ion and the p-wave states 
of the incident ele ctron. The coup ling is responsible for the fast 
DR rate (Kokoouline et al.ll2001h in H^. In the present model, 
only ffie p-wave electronic states are included because other par- 
tial waves have a much smaller effect on the DR probability : the 
s-wave states have no E-type character and, therefore, are only 
weakly coupled to the dissociative electronic states of H3; d- 
wave electronic states are coupled to the Zi-vibrational modes, 
but the coupling is raffier small because the cZ-wave of the in- 
cident electron does not penetrate sufficiently close to the ionic 
core owing to the cZ-wave centrifugal potential barrier. 

All three internal vibrational coordinates are taken into ac- 
count. Vibrational dynamics of the ionic core are described using 
the hyper- spherical coordinates, which represent the three vibra- 
tional degrees of freedom by a hyperradius and two hyperangles. 
The hyperradius is treated as a dissociation coordinate that rep- 
resents uniformly the two possible DR channels, the three-body 
(such as H-hH-i-H) and two-body (such as H2-1-H). Although the 
initial vibrational state of the ion is the ground state, after re- 
combination with the electron, other vibrational states of the 
ionic target molecule can be populated. Therefore, in general, 
many vibrational states have to be included in the treatment. In 
particular, the states of the vibrational continuum have to be in- 
cluded, because only such states can lead to the dissociation of 
the neutral molecule. The vibrational states of the continuum are 
obtained using a complex absorbing potential placed at a large 
hyperradius to absorb the flux of the outgoing dissociative wave. 

Since the rovibrational symmetry is D3/, for and D^ and 
C21. for H2D^ and D2H^, the rovibrational functions are clas- 
sified according to the irreducible representations of the corre- 
sponding symmetry groups, i.e. A'p A", Aj, Aj, E', and E" for 



D^h and Al, A2, Bi, and B2 for C21.. We use the rigid rotor approx- 
imation, i.e. the vibrational and rotational parts of the total wave 
function are calculated independently by diagonalizing the cor- 
responding Hamiltonians. In our approach, the rotational wave 
functions must be obtained separately for the ions and for the 
neutral molecules. They are constructed in a different way for 
the D3/, and C^v cases. The rotational eigenstates and eigenen- 
ergies of the ^3;, molecules are symmet ric top wave functions 
(see, for example [Bunker & Jenseniri998h . They can be obtained 
analytically if the rotational constants are known. The rotational 
constants are obtained numerically from vibrational wave func- 
tions, i.e. they are calculated separately for each vibrational 
level of the target molecule. The rotational functions for the C2v 
ions are obtained numerically by diagonalizing ffie asymmetric 
top Hamiltonian dBunker & Jensenlll998llKokoouhne & Greend 
2005). 

Once the rovibrational wave functions are calculated, we 
construct the electron-ion scattering matrix (5 -matrix). The 
5 -matrix is calculated in the framew ork of quantum de- 
fect theory (QDT) (see, for example, 'Chang & Fand 1 19721 
Kokoouline & Greene 2003b, 2005) using the quantu m defect 
parameters obtained from the ab initio calculation dMistrik et al.l 
2001). The constructed scattering matrix accounts for the Jahn- 
Teller effect and is diagonal with respect to the different irre- 
ducible representations F and the total angular momentum of 
the neutral molecule. Thus, the actual calculations are made sep- 
arately for each F and A^. Elements of the matrix describe the 
scattering amplitudes for the change of the rovibrational state 
of the ion after a collision with the electron. However, the S- 
matrix is not unitary due to the presence of the dissociative vi- 
brational channels (i.e. continuum vibrational states of the ion, 
discussed above), which are not explicitly listed in the computed 
5 -matrix. The "defect" from unitarity of each column of this 
5 -matrix is associated with the dissociation probability of the 
neutral molecule formed during the scattering process. The dis- 
sociation probability per collision is then used to calculate the 
DR cross-sections and rate coefficients. 

The nuclear spin states are characterized by one of the A 1 , 
A2, or E irreducible representations of the symmetry group S 3 
for D3;, molecules and by the A or B irreducible representations 
of the symmetry group 5 2 for C2,, molecules. The irreducible 
representation F„5 of a particular nuclear spin state determines 
its statistical weight and is related to the total nuclear spin / of 
the state. Here, / is the vector sum of spins i of identical nuclei. 

For H3 , the F„j = A[ states {A'^ and A'^ rovibrational states) 
correspond to / = 3/2 (ortho); the F„s = E states (E' and E" 
rovibrational states) correspond to 1-1/2 (para). The statistical 
ortho:para weights are 2 : 1 . 

For H2D^, the r„j = A states (Bi and B2 rovibrational states) 
correspond to / = 1 (ortho); the F„j = B states (Ai and A2 rovi- 
brational states) correspond to / = (para). The statistical or- 
tho :para weights are 3:1. 

For D2H^, the r„j = A states (Ai and A2 rovibrational states) 
correspond to / - 0,2 (ortho); the r„j = B states (Bi and B2 
rovibrational states) correspond to / = 1 (para). The statistical 
ortho:para weights are 2 : 1 . 

Finally, for D3 , the r„ s = Ai states (A'j and A" rovibrational 
states) correspond to / = 1,3 (ortho); the r„j = A2 states (Aj and 
Aj rovibrational states) correspond to / = (para); the r„5 = E 
states (£" and E" rovibrational states) correspond to / = 1,2 
(meta). The statistical ortho:para:meta weights are 10:1:8. 

Figures dB.lt . (|B.2t . dB.3t . and (|B.4t summarize the ob- 
tained DR thermal rate coefficients calculated separately for each 
nuclear spin species of the four H3 isotopologues and the numer- 
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Fig.B.l. Theoretical DR rate coefficients as functions of tem- 
perature for the ortho- and para-species of . The figure also 
shows the species-averaged rate coefficient. For comparison, we 
show the rate coe fficient obtained from measurements in the 
TSR storage ring bv Kreck el et alj (l2005l) and the analytical de- 
pendence for the coefficient used in earlier models of prestellar 
cores by FPdFW. 
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Fig. B.2. Theoretical DR rate coefficients for the ortho- and 
para-species of HaD^, as functions of temperature. The rate co- 
efficient averaged over the two species and the analytical expres- 
sion used in earlier models are also shown. 



ical values are listed in Table IB. 11 For comparison, the figures 
also show the analytical dependences used in previous models 
of prestellar core chemistry (FPdFW). As one can see, the rates 
for diff'erent nuclear spin species are similar to each other (for 
a given isotopologue) at high temperatures. However, for lower 
temperatures, the rates for diff'erent ortho/para/meta-nuclear spin 
species significantly differ from each other. The diff'erence in be- 
havior at low temperatures is explained by different energies of 
Rydberg resonances present in DR cross-sections at low electron 
energies. The actual energies of such resonances are important 




Fig. B.3. Theoretical DR rate coefficients for the ortho- and 
para-species of D2H^ as a function of temperature. The rate co- 
efficient averaged over the two species and the analytical expres- 
sion used in earlier models are also shown. 
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Fig. B.4. Theoretical DR rate coefficients for the ortho-, para-, 
and meta-species of D^. The rate coefficient averaged over the 
two species and the one used in earlier models are also shown. 



for the thermal average at temperatures below or similar to the 
energy difference between ground rotational levels of different 
nuclear spin species. At higher temperatures, the exact energy of 
the resonances is not important. The averaged rate is determined 
by the density and the widths of the resonances, which are sim- 
ilar for all nuclear spin species over a large range of collision 
energies. 
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.65 


8, 


.57e- 


-09 


7, 


.66e- 


■09 


1, 


.45e- 


-08 


1, 


.15e- 


-08 


1, 


.Ole- 


■08 


8, 


.66e- 


-09 


2, 


.71e- 


■09 


3, 


.08e- 


■09 


2, 


.90e- 


-09 


4310, 


.04 


7, 


.75e- 


-09 


6, 


.91e- 


■09 


1, 


.30e- 


-08 


1 


.03e^ 


-08 


9, 


.04e- 


■09 


7, 


.76e- 


-09 


2, 


.40e- 


■09 


2, 


.72e- 


■09 


2, 


.57e- 


■09 


4758 


.74 


6 


.99e- 


-09 


6, 


.23e- 


■09 


1 


.16e^ 


-08 


9 


.26e^ 


-09 


8, 


.lOe- 


■09 


6 


.95e^ 


-09 


2, 


.13e- 


■09 


2, 


.41e- 


■09 


2, 


.27e- 


■09 


5254 


.15 


6 


.31e- 


-09 


5, 


.61e- 


■09 


1 


.04e^ 


-08 


8 


.31e^ 


-09 


7, 


.26e- 


■09 


6 


.23e^ 


-09 


1 


.89e- 


■09 


2, 


.14e- 


■09 


2, 


.02e- 


■09 


5801 


.13 


5, 


.70e- 


-09 


5, 


.06e- 


■09 


9 


.39e- 


-09 


7, 


.47e- 


-09 


6, 


.51e- 


■09 


5 


.58e- 


-09 


1 


.68e- 


■09 


1 


.90e- 


■09 


1, 


.79e- 


■09 


6405, 


.05 


5, 


.15e- 


-09 


4, 


.58e- 


■09 


8, 


.45e- 


-09 


6 


.73e- 


-09 


5, 


.85e- 


■09 


5, 


.02e- 


-09 


1, 


.50e- 


■09 


1, 


.69e- 


■09 


1, 


.60e- 


■09 


7071, 


.85 


4, 


.67e- 


-09 


4, 


.14e- 


■09 


7, 


.63e- 


-09 


6 


.08e^ 


-09 


5, 


.27e- 


■09 


4, 


.52e- 


-09 


1, 


.34e- 


■09 


1, 


.51e- 


■09 


1, 


.43e- 


■09 


7808 


.06 


4, 


.24e- 


-09 


3, 


.77e- 


■09 


6 


.91e^ 


-09 


5 


.50e^ 


-09 


4, 


.77e- 


■09 


4 


.08e^ 


-09 


1 


.20e- 


■09 


1 


.36e- 


■09 


1 


.29e- 


■09 


8620 


.92 


3, 


.87e- 


-09 


3, 


.43e- 


■09 


6 


.28e^ 


-09 


5 


.Ole^ 


-09 


4, 


.33e- 


■09 


3, 


.71e^ 


-09 


1 


.09e- 


■09 


1 


.23e- 


■09 


1 


.16e- 


■09 


9518 


.40 


3, 


.55e- 


-09 


3, 


.14e- 


■09 


5, 


.73e- 


-09 


4, 


.57e- 


-09 


3, 


.95e- 


■09 


3, 


.38e- 


-09 


9 


.86e- 


■10 


1, 


.lle- 


■09 


1, 


.05e- 


■09 


10509 


.31 


3, 


.26e- 


-09 


2, 


.89e- 


■09 


5, 


.26e- 


-09 


4, 


.20e- 


-09 


3, 


.61e- 


■09 


3, 


.09e- 


-09 


8, 


.99e- 


■10 


1, 


.Ole- 


■09 


9 


.60e- 


■10 


11603 


.38 


3, 


.Ole- 


-09 


2, 


.67e- 


■09 


4, 


.85e^ 


-09 


3, 


.87e^ 


-09 


3, 


.33e- 


■09 


2, 


.85e- 


-09 


8, 


.25e- 


■10 


9 


.29e- 


■10 


8, 


.80e- 


■10 


12811 


.35 


2, 


.80e- 


-09 


2, 


.48e- 


■09 


4 


.50e^ 


-09 


3 


.58e^ 


-09 


3, 


.08e- 


■09 


2, 


.63e^ 


-09 


7, 


.61e- 


■10 


8 


.56e- 


■10 


8 


.12e- 


■10 


14145 


.07 


2, 


.61e- 


-09 


2, 


.31e- 


■09 


4 


.19e^ 


-09 


3, 


.34e^ 


-09 


2, 


.86e- 


■09 


2, 


.45e^ 


-09 


7, 


.05e- 


■10 


7, 


.94e- 


■10 


7, 


.53e- 


■10 


15617 


.64 


2, 


.45e- 


-09 


2, 


.16e- 


■09 


3, 


.92e- 


-09 


3, 


.12e- 


-09 


2, 


.67e- 


■09 


2, 


.29e- 


-09 


6 


.58e- 


■10 


7, 


.40e- 


■10 


7, 


.Ole- 


■10 


17243, 


.51 


2, 


.30e- 


-09 


2, 


.04e- 


■09 


3, 


.68e- 


-09 


2, 


.94e- 


-09 


2, 


.51e- 


■09 


2, 


.14e- 


-09 


6 


.16e- 


■10 


6 


.93e- 


■10 


6 


.57e- 


■10 


19038 


.65 


2, 


.18e- 


-09 


1, 


.92e- 


■09 


3, 


.48e- 


-09 


2, 


.77e- 


-09 


2, 


.37e- 


■09 


2, 


.02e- 


-09 


5, 


.80e- 


■10 


6 


.52e- 


■10 


6 


.19e- 


■10 


21020 


.66 


2, 


.07e- 


-09 


1, 


.83e- 


■09 


3 


.30e^ 


-09 


2, 


.63e^ 


-09 


2, 


.24e- 


■09 


1 


.92e^ 


-09 


5 


.49e- 


■10 


6 


.17e- 


■10 


5, 


.85e- 


■10 


23209 


.01 


1 


.97e- 


-09 


1, 


.74e- 


■09 


3, 


.15e^ 


-09 


2, 


.51e^ 


-09 


2, 


.13e- 


■09 


1 


.82e^ 


-09 


5 


.22e- 


■10 


5 


.86e- 


■10 


5 


.56e- 


■10 



Table B.l. Dissociative recombination rates of , H2D*, D2H*, and for each 
individual nuclear spin state species. 



